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ABSTRACT 

We use a sample of blue horizontal branch (BHB) stars from the Sloan Digital Sky Survey Data Release 7 to 
explore the structure of the tidal tails from the Sagittarius Dwarf Galaxy. We use a method yielding BHB star 
candidates with up to ^ 70% purity from photometry alone. The resulting sample has a distance precision of 
roughly 5% and can probe distances in excess of 100 kpc. Using this sample, we identify a possible extension to 
the trailing arm at distances of 60-80 kpc from the Sun with an estimated significance of at least 3.8(t. Current 
models predict that a distant 'returning' segment of the debris stream should exist, but place it substantially 
closer to the Sun where no debris is observed in our data. Exploiting the distance precision of our tracers, we 
estimate the mean line-of-sight thickness of the leading arm to be ~3 kpc, and show that the two 'bifurcated' 
branches of the debris stream differ by only 1-2 kpc in distance. With a spectroscopic very pure BHB star 
subsample, we estimate the velocity dispersion in the leading arm, 37 km s"' , which is in reasonable agreement 
with models of Sgr disruption. We finally present a sample of high-probability Sgr BHB stars in the leading 
arm of Sgr, selected to have distances and velocities consistent with Sgr membership, to allow further study. 
Subject headings: Galaxy: halo — stellar content — galaxies: dwarf — interactions 



1. INTRODUCTION 

Tidal debris from dwarf galaxies and stellar clusters 
dissolving in the Milky Way potential are an impor- 
tant contributor to the stellar halo of the Milky Way 
(e.g., fSearie & Zinn 1978; Ibata et al. 1994; Bullock et^ 
2001; Bullock & JohnstonI 120051: IBelokurov et al.i 1200631 
Bell et al. 2008i ). In recent years, many elongated sub- 
structures have bee n found i n the s tellar halo of the 
Milky Way (e.g., 'Ibata et al.' l99f, "^OOS"; 'Yannv et alj 
2003; Grillmair & Johnson 2006; Grillmair & Dionatos 20061 
Grilknaiil |2006; Belokurov et al. 2007) and around other 
nearby galaxies suc h as Andromeda (e.g., Ibata^et al. 2001al 
iMcConnachie et all 12009). an d a number of e xternal galax- 
ies (e.g., N GC 891; Mouhcine et al.' "201^, NGC 5907; 
IZheng et al . 1999; Martmez-Delgado et al. 2008, 2010) show- 
ing that the build-up of stellar halos through accretion of satel- 
lite galaxies is a common phenomenon. Besides the general 
implications such stellar satellite debris has for building and 
testing the galaxy formation paradigm, the detailed investiga- 
tion of the individual structures provides important informa- 
tion about the specific formation history of individual galax- 
ies. The spatial distribution and kinematics of the tidal de- 
bris of dwarf galaxies or globular clusters is also an impor- 
tant source of informati on about the g ravitat ional p otential 
of the Milky W a y (e.g.. Fjohnston et alj[T999l iHelmiilllOOia : 
Law et al.' '2005^. 'Fellhauer et all 120061: iKoposov et alJ 12009 : 
Law & Majewski 2010a; Pena rrubia et al.ll2010al) ! 

In this context, the Sagittarius stellar stream (Sgr), the most 
massive stellar stream arou nd the Milk y Way, is a central case 
study. Discovered in 1994 (ilbata et al.iri994f) . the tidal tail has 
been charted acro ss more than one full w rap arou nd the Milky 
Way in M-giants (iMajewski et alJlToOSL see alsolYanny et alJ 
l2009l) . main sequence stars teelokurov et al.l 1200631) . clus- 



'Monaco et al." '2003; Clewlev & J3rvis| 120061: lYannv et all 
2009; Niederste-Ostholt et al. 2010). The sp3ti3l tightness of 
the stream in combination with its full 360° span m3kes it 3n 
important probe of the potential (e.g.,'Helmi & Whitdll99^ 
Moore etal. 1999; Ibata etal. 2001c, 2002; Johnston et 3 
2002^2005; Helmi 2004a; Lewis & Ibata 2005; Binnev 20oC 
of the disruption process (Ibata et al. 2001b; Helmi & Whi^ 
1200 U [Penarrubi a et 3LI 120 10b). and of the impact of pop- 
ulation gradients 3nd cluster contents of the Sgr dwarf 
on th e properties of the tai l (e.g.. iDa C osta & Armandrof) 



1995*; 'Maiewski et 3I 



120031: lM3rtmez-Delg3do et 3L ,200- 
L3W & M3iewskill20I0bl) . 



ters (e.g. |Bellazzini et al. 2003] and re ferences therein), 3nd 
blue horizontal branch (BHB) stars (iNewberg et 3I.I l2003l 



B_ellazzini et al. 2003: 

Despite the wealth of observational d3t3, models of the 
stream have failed so far to match all the observational con- 
straints by quite a margin. To explain the observations 
different gakxy p otentisls h3ve been invoked, with argu- 
ments for prolate (iHelmil l2004bl: iLaw et 3l.ll2()05h. sp heric3l 
( lFellh3 ueret al.' '2006*). oblate ([Johnston et all 120051) or tri- 
axi3l (fLaw & Maiewski 20103*) dark matter potentials. To 
explain some striking fe3tures, s uch 3S the 'bifurc3tion' 
(Belokurov et 3I. 20063), Pen3rr ubi3 et 3I.I (f20 10b) invoked 
th3t the progenitor of the Sgr stre3m m3y h3ve been 3 rot3ting 
disk g3l3xy r3ther th3n 3 pressure-supported dw3rf g3l3xy 3S 
3ssumed by most previous models. However, no single mod- 
els seems to expl3in 3II p3rts of the stre3m while it is 3lso not 
entirely cle3r th3t 3II the overdensities found in the pkne of 
the Sgr stre3m 3re 3ctu3lly remn3nts of the S3me progenitor 
A more precise 3nd more complete empiric3l picture of the 
Sgr stre3m could be cruci3l in ckrifying this issue, 3nd this 
constitutes the centr3l go3l of the present p3per 

In recent studies of the Sgr stre3m, there h3s been incre3sed 
3ttention tow3rd BHB st3rs 3S 3 tr3cerpopul3tion. Due to their 
rel3tive brightness they C3n be observed out to ~' 100 kpc in 
the stell3r h3lo of the Milky W3y using Slo3n Digit3l Sky Sur- 
vey (SDSS) d3t3. However, to t3ke full 3dv3nt3ge of 3re3 cov- 
er3ge of surveys such 3S SDSS, the identific3tion of these st3rs 
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needs to be done with photometric data al one. Many pubhca- 
tions based their selection on color b oxes (lYannv et al.ll2000l 
120091: iNiederste-Ostholt et al.1 120101 ) that included a signifi- 
cant contamination from other blue stars (primarily blue strag- 
gler (BS) stars). Such contaminants can dominate in number, 
and are 1-2 mag fainter in absolute magnitude, confusing the 
interpretation of halo structure using such samples. 

In this paper, we use SDSS data in the North Galactic Cap 
to study Sgr tidal debris. We choose color-selected BHB 
star candidates as sparse tracers of the ancient, metal poor 
populations with well-defined absolute magnitudes, that are 
^3-4 magnitudes brighter than the densely populated main- 
sequence tu rn-off (MSTO) stars. Going beyond other recent 
studies (e.g.. lYannv et al.l2009t:lNiederste-Ostholt et al.l20"Tol) 
of the Sgr system in BHB stars we use a refined selection 
technique based on a spectroscopic training sample which re- 
duces the contamination by other stellar populations (Bell et 
al. 2010). We show empirically that the distance uncertain- 
ties in our sample are small, of the order of 5%. We use 
these stars to chart out the Sgr stream, focusing on three is- 
sues: delineating the distant (> 50 kpc) overdensities that 
may be associated with the Sgr trailing arm, on constraining 
and measuring the thickness of the leading arm, and on pre- 
senting a sample of high-probability Sgr BHB star candidates 
with positions and velocities consistent with Sgr membership 
for further study. Furt hermore we are explore the bifurcation 
that has been found bv lBelokurov et all ( 12006a ) perpendicular 
to the orbital plane of the stream and its appearance in BHB 
stars. 

2. DATA 

2.1. Blue Horizontal Branch Stars 

For this study, we use Data Release 7 ('Abazaji an et al.l 
12009) of the SDSS to probe the Sagittarius stellar stream with 
BHB stars. The SDSS is an imaging and spectroscopic survey 
that has so far mapped a little over ~ 1 /4 of the sky. Imaging 
data are produced simultan eously in five p hotometric bands, 
namely, u, r, and z (Fukugita et al.' '1 9961 : iGunn et al] 
[1998; Hog g et alJ 120011 Gunnetal. 2006). The data are 
processed through pipelines to measure photometric and 
astrometric properties (iLupton et a 1[T999[ IStoughton etal J 
2002; Smi th et al.ll2002t iPier et al.ll2003t llvezic e t al. 2004: 
Tucker et al T and to se lect targets for spe ctroscopic 

follow-up iBlanton et al. 2003: IStrauss et al.ll2002l) . 

The horizontal branch is populated by stars which have 
developed past the main sequence stage and are now burn- 
ing helium in their cores and hydrogen in the shell. BHB 
stars have the dual advantages of a high luminosity (al- 
lowing probing of the Milky Way halo to > 100 kpc), and 
have a small intrinsic spread in absolute magnitudes. Their 
main disadvantage is that the selection of a clean sam- 
ple of BHB stars is challenging from photometry alone. 
While broad cuts in u- g and g - r are sufficient to iso- 
late BHB stars and other A-type stars (expecte d to be BS 
stars; iPreston & SnedenI 120001 ISirko et al.l [2004b from low- 
redshift quasars and white dwarfs, distinguishing BHB stars 
from the BS contaminants is considerably more challenging 
(e.g., Kinman et al. 1994; Wilhelm et al. 1999; Clewlev et aU 
2002; Sirko et al. 2004; Kinman et al. 2007; Xue et al. 2008; 
Smith et al..2010.) . Previous works have us ed broad color cuts 
des igned to mitigate this c ontamination (lYannv et al.l 120091 
andlNiederste-Ostholt et al.''2Q10 used the selection in Figure 
10 of [Yannv et al...2000; Sirko et al.. 2004 used a different cut 
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Fig. 1 . — Photometric BHB stai' selection. All ca ndidate BHB star s shown 
have with spectroscopic classification (following IXue et al.l 12008). Stars 
spectrally classified as BHB stars are shown in blue; other stars are shown 
in red. Contours show regions in (u—g) vs. (g—r) space where the fraction of 
targets classified as BHB star's exceeds 10%, 30% and 50%, respectively; the 
50% contour is thicker than the others. In what follows we focus primarily 
on BHB star candidates whose colors fall within the 50% contour. 

for their faint sample of BHB candidates). Yet, these methods 
all suffer from very substantial contamination from BS stars. 

Spectroscopy permits a fairly clean separation of BHB 
stars from BS stars o n the basis of s urface gravity dependent 
B aime r line profiles. IXue et all (12008 ), following ISirko et alj 
( 120041) . use a two-stage cut to distinguish BHB from BS 
stars. First, stars in the color box 0.8<m-^<1.6 and 
-0.5 < g-r < 0.0 with a relatively low line width and low 
flux in the line core relative to the continuum are chosen (this 
reduces contamination to ^ 50%). Then, a Sersic profile is fit- 
ted to the Balmer lines. By combination of these two criteria, 
a > 90% pure sample of BHB stars is isolated. Unfortunately, 
SDSS spectroscopy of BHB stars (mostly from SEGUE) is 
limited to certain areas of sky, and only relatively bright BHB 
stars are targeted, meaning that BHB stars more distant than 
> 50 kpc are not well-probed by the SDSS. 

Therefore, we have re-addressed the issue of photometric 
selection of BHB star candidates (described in full in Bell et 
al. 2010). We use the spectroscopic classifications of < 18 
stars from Xue et al. (2008) as a training set. We calculate the 
probability of a star in the color box 0.8 < m-^ < 1.6 and 
-0.5 < g- r < 0.0 being a BHB star from this training set 
(Figure[T|l. Blue data points show stars that are very likely to 
be BHB stars on the basis of their spectra (a contaminatio n of 
much less than 1 0% has been argued by IXue et al.ll2008l and 
ISirko et ani2004l for < 18). The thick contour outlines the 
region of color-color space where the fraction of BHB candi- 
dates that are spectroscopically-classified BHB stars is > 50% 
and there were more than 16 stars in a bin of .025 x 0.04 mag. 
Applying this selection to the SDSS DR7 ( lAbazaiian et alj 
"2009), we obtain a candidate sample with 389,785 stars within 
the 0.8 < M-g < 1.6, -0.5 < g-r < 0.0 color box. In the 



The Structure of Sgr in BHB Stars 



3 



following we apply a lower probability limit of 50% for the 
photometric sample reducing the sample size to 28,270 stars. 
Tests show that this '> 50%' probability sample isolates half 
of the mg < 18 BHB star population, with a contamination 
of 20-30%. Performance at fainter limits is expected to de- 
grade gradually, with increasing incompleteness and contam- 
ination (at irig ^ 20 roughly 1/4 of BHB stars are expected 
to be kept, and contamination may be as severe as 50%; Bell 
et al. 2010). We will later test the influence of changing the 
probabilit y cu ts (and therefore completeness/contamination) 
in Section im 

2.1.1. Kinematic Sample 

The radial velocity sample, which is a sub-sample of the 
photometric sample, was selected based on the spectra as de- 
scribed above offering a much higher BHB purity > 90% than 
the method applied on the stars with photometry only. To not 
unnecessarily restrict the sample size, we use the full radial 
velocity sample in these cases and ignore for these stars the 
probabilities which were assigned based on their colors (i.e., 
we do not use the lower probability limit of 50% mentioned 
above). The total sample size is 5233 stars, of which 807 are 
located in the Sgr plane (see Section |2.1.2| |. From these 807 
stars 616 would fulfill the 50% probability criterion, giving 
a success rate of spectroscopic BHB stars in this selection of 
76%. Throughout this paper the radial velocities are given in 
the Galactic standard of rest, which are the heliocentric ra- 
dial velocities corrected for the Galactic rotation assuming a 
rotation velocity of 220 km s"' for the local standard of rest 
and (H-10.0,H-5.2,H-7.2)kms"' for the solar motion where the 
directions are defined as pointing towards the Galactic center, 
in t he direction of ro tation and towards the north Galactic Pole 
(see lXue et alll2008l for details). 

2.1.2. Sagittarius in a Galactic Plane 

For much of our analysis, we focus on stars in the presumed 
orbital plane of the Sgr stream only. We define this 'Sagit- 
tarius plane' to encompass the Sgr stream and the Galactic 
Center; this is presumably close to the orbital plane of the Sgr 
stream. To ensure consistency with models, we use the same 
pole as the Two Micro n All Sky Survey (2MASS) papers (e.g. 
lMaiewskietaIll2003h at {l,b) = (273°.8 -13°.5). Stars are 
considered to be in the plane if they lie within ± ^ 27° of 
this plane; this definition naturally yields not a plane but a 
wedge, whose physical thickness increases with distance from 
the Sun. Stars are projected onto this plane by conserving the 
distance to the Sun (i.e., the plane is a projection of shell seg- 
ments onto the plane). The Sagittarius plane defined here in- 
cludes 73,066 stars (6905 with a BHB star probability greater 
than 50%) from the total 389,785 stars (28,270 with a BHB 
star probability greater than 50%) in the SDSS volume that 
are inside the color box. 

2.1.3. Empirical Distance Uncertainties 

As distance precision for the BHBs plays an important role 
for our analysis, we use several known globular clusters and 
dwarf spheroidals to determine both the statistical and sys- 
tematic uncertainties of the distance determinatioifl A color- 
magnitude diagram (CMD) for one of the dwarf spheroidals 

The distances were derived using the g — r-dependent Mg calibration 
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Fig. 2. — BHB section of the color magnitude diagram for the Ursa Minor 
dwarf spheroidal galaxy. The diamonds represent the ugr-color selected BHB 
star candidates with > 50% BHB probability; the rest of the sample with 
lower probabilities is shown as crosses. It can be clearly seen that there is 
a trend towards fainter magnitudes for bluer colors; the g-r dependent Mg 
cahbration of Sirlco et al., (2004) follows this trend closely. The vertical line 
shows the position of the color cut applied to distinguish between red and 
blue BHB stars. 

is shown in Figure |2] The typical shape of the blue horizon- 
tal branch shows a nearly horizontal part at redder colors and 
a gradual trend towards fainter magnitudes at the blue end, 
as can be seen in Figure |2] Overall this trend causes an in- 
crease in the magnitude spread and therefore distance mea- 
surement unce rtainties toward bluer colors. As we will show 
in Section ITT] the Sagittarius stream shows a larger concen- 
tration of 'red' BHB stars. This indicates that studying the red 
BHB stars separately can have two benefits compared to look- 
ing only at the sample as a whole. /) The signal strength for 
the stream will increase, and //) the uncertainties introduced 
by the deviations from the horizontal shape of the horizontal 
branch can be reduced. 

Therefore, we divide the sample into a blue and a red part 
for further analysis. The g-r value at which we apply the 
cut throughout this paper is illustrated in Figure |2]by the ver- 
tical line. This cut is chosen to divide the bright stars of the 
sample (^ < 18.5 mag) in equally populated halves. We de- 
termine the statistical error of the distance measurement for 
BHB stars by measuring the spread of their distance moduli 
within one cluster (whose line-of-sight extent is negligible). 
The distance modulus distribution for the objects is shown in 
Figure |3] We fit Gaussians to these distributions and use the 
standard deviation for estimating the statistical distance un- 
certainty 6D/D. We measure the mean value and the standard 
deviation for both the red part and the > 50% BHB probabil- 
ity sample (see TablefTli. The distribution in distance modulus 
of red and blue stars is also shown in Figure|3]indicated by the 
blue and red shaded areas. The results are shown in Table [T] 
The mean statistical distance uncertainty for the objects listed 
here is 4% for the > 50% sample and 6% for the full sample. 

Comparison with prior distance determinations (see Table 
[U showed a systematic underestimation of the distances in 
our results. This effect is of the order of 4% in distance, but 
also includes some variance which is probably also partly due 
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Fig. 3. — Distance accuracy for individual BHB stars. Tiie panels show the distance modulus distribution of photometiically selected SDSS BHB stars in four 
clusters and dwarf spheroidal galaxies. For each object, a Gaussian is fitted to the full sample (unfilled histogram) to measure the width of the distiibution as an 
estimate of the distance uncertainty. The gray histogram shows only the stars with a BHB probability greater than 50% (Figure[T], the ones filled with vertical red 
and horizontal blue lines give the distribution for the red and the blue pait of the hoiizontal branch, respectively. Cle^ly these clusters have a larger concentration 
of red BHB stars. Results of the fit to both the > 50% BHB probability sample and the red part are shown in Table[T] The statistical distance uncertainty SD/D 
resulting from the fits to the full sample is also given in the plots, being a bit higher for the full sample and a bit lower for the > 50% sample than the expected 
value of SD/D ~ 0.05 tSirko et al.l2004l) . 



to the fact that the literature values were determined with dif- 
ferent methods. 

With this test we cannot probe uncertainties in the distance 
determination that arise from a spread i n metalhcity . The 
metallicity-dependentBHB star models of lDotter et alJ (120071) 
and Dotteret al. (2008) indicate a significant contribution to 
the distance uncertainties by a range of metallicities in the 
halo BHB stars. The overall uncertainty accounting for a com- 
bination of the scatter we see in single metallicity populations 
and the contribution of a scatter introduced by having a va- 
riety of metallicities is estimated to be less than 10% in Bell 
et al. (2010). In what follows we account only for the un- 
certainty which was estimated using single metallicity popu- 
lations, which may underestimate the overall distance uncer- 



tainties (5% vs. < 10%). 

As a comparison da ta set we use M giants from the 2MASS 
dSkrutskie et al.l 120061) to compare the distance scale of our 
BHB star data set in relation to other stellar populations, 
which were used for studying the Sgr stellar stream. In par- 
ticular, this M giant data set was also used as the basis for 
the models we will compare to later M giants can be used as 
distance indicators out to large distances making them a good 
stellar population for studying the Sgr system (especially in 
the near infrared). Due to the complete coverage of the sky, 
it is possible to observe the stellar stream along its whole or- 
bital path. A disadvantage of M giants as distance indicators is 
their rather large distance uncertainties (argued to be ^ 17%; 
iLaw et al.1 120051) . and a likely distance offset with the BHB 
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TABLE 1 

Distances to Clusters and Dwarf Spheroidals 



Object 


{m-M)prob>0.i 


O" pmb>O.S 


{m-M)red 




Literature Values 


SD/D 


Rel. Dist. Offset 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 






NGC 5024 


16.26 


0.09 


16.27 


0.14 


16.31" 


0.04 


0.02 


Bootes 


18.95 


0.08 


18.95 


0.12 


18.94 ±0.14'' 


0.04 


0.01 


Ursa Minor 


19.17 


0.07 


19.15 


0.07 


19.32 ±0.12*^ 


0.03 


0.07 


Sextans 


19.56 


0.12 


19.57 


0.11 


19.75 ±0.13'' 


0.05 


0.06 



Note. — Mean and standard deviation in distance modulus for four clusters and dwarf galaxies given for a subsample with a BHB star probability > 50% 
and for the red BHB stars only. The second to last column gives the inferred distance uncertainty SD/D for the > 50% sample of the color-selected BHB stars. 
The mean value is 0.04, which is a bit lower than the mean value of the full sample of 0.06. As hterature values we give the mean values and standard deviations 
from a number of studies as listed. In the last column, the relative distance offset is given which was calculated using the literature values listed in the table. 

jHarrisll tT^ . 

^elokurov et al. (2006b); Dall'Ora et al. (2006); Sieael (2006); de Jona et al. (2008). 
Wighell & B urke 1 1999J; Bellazzini et al. ( 2002); Carrera et ah (2002.) ; .Tammann et al] fSOOg) . 
HMateo et all {l99^ ; lLe"e et alj <2003l) ; ITammann etaH f200a) . 



-100 



-50 - 



- 



50 - 



100 




-50 
distance from Sun 



n kpc 



Fig. 4. — Distiibution of M giants from 2MASS in the Sagittarius orbital 
plane is shown as r ed do ts with a gray BHB star probability map in the back- 
ground (see Section[3T). For both the BHB stars and the M giants we exclude 
the inner 20 kpc. The M giant population in the leading arm (-1-20, -h35) ap- 
pears closer to the Sun than the I3HB stars, which we inteipret as an 8% 
mismatch in the distance scale between the two different populations. Fur- 
thermore, the leading arm appears to be much more compact in width in BHB 
stars than in M giants, presumably a reflection of the substantially larger dis- 
tance uncertainties in the M giants. 

and literature distance scales. 

We derive a sample of M giants fro m the full 2M ASS cata- 
log following the method described in lMajewski et al. (2003) 
for which we show the distribution in the plane of the Sgr 
stellar stream in Figure |4] The comparison with the BHB star 
population also shown in this plot reveals a distance offset 
between the two populations with the M giants being about 
~ 8% closer to the Sun than the BHB stars in the leading arm 
region of Sgr As the distances of M giants are less well de- 
termined than those of BHB stars we also see a difference in 
the width of the leading arm in the different populations; the 
width seen in BHB stars is only ^ 40% of that seen in M gi- 
ants as it appears in the samples presented here. Obviously 
this mismatch will propagate through to the models based on 
M giant observations (e.g. . .Law et al.,.2005:.Law & Maiewski 



l2010al) . so that we are expecting to see this mismatch to some 
degree in the comparison to these models. Note that we do 
not adjust our distance scale (or those of other data or mod- 
els) to account for possible distance offsets in either case (the 
^ 4% mismatch between the BHB distance scale and the lit- 
erature determinations, or the ~ 8% mismatch between the 
BHB and M giant distance scales). In particular, this means 
that throughout the paper different data sets or models shown 
in the same plot can have different distance scales. Note that 
the offset to the distance scales of the two comparison sam- 
ples have opposite directions, while the clusters and dwarf 
spheroidals have systematically larger distances, the M gi- 
ants have smaller distances compared to our BHB star sam- 
ple. This implies an even larger offset between these distance 
scales of about 12%. Evidently a better characterization of 
the M giant distance scale would be of importance for a di- 
rect comparability of different stellar populations as distance 
indicators. 

2.2. N-Body Models for the Sgr Stream 

We will compare our BHB maps with simulations of the 
evolution of the Sagittariu s dwa r f spheroidal in the Milky 
Way potential (Law ej a]] 120051; iLaw & Maiewskil 1201 Oat 
iPefiarrubia et alJ l2010bl) and summarize these models here. 
The lLaw et al.l (l2005h models adopt a smooth, rigid potential 
representing the Milky Way, which consists of a Miyamoto- 
Nagai disk, a Hemquist spheroid, and an axisymmetric log- 
arithmic halo of different flattenings: q = 0.9 (oblate), 1.0 
(spheroidal ), and 1.25 (pro l ate). W e will also use a new 
model by Law & Maiewskil ( 12010a ) for comparison, which 
is based on a triaxial dark matter halo with a minor/major 
axis ratio (c/a)$ = 0.72 and a intermediate/major axis ratio 
(b/a)<j, = 0.99 at radii 20 < r < 60 kpc. This corresponds to a 
nearly-oblate ellipsoid whose minor axis is contained within 
the Galactic disk plane and approximately aligned with the 
line of sight to the Galactic Center In both model genera- 
tions, the Sagittarius dwarf itself is represented by lO*" self- 
gravitating particles. All of the models were constructed to fit 
the system of the Sagittarius stellar stream as seen in 2MASS 
M giants. To account for the photometric distance errors of 
the M giant sample, a artificial random distance error of 17% 
was applied to the simulated debris particles. Followin g the 
suggestion of a triaxial halo, IPefiarrubia et al.l (1201 Obi pre- 
sented a model which does not assume a pressure-supported 
dwarf spheroidal galaxy as the progenitor of the Sgr stellar 
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Fig . 5 . — Density maps of BHB stars in a plane which includes the Sgr t idal 
tail and the Galactic Center. The maps ai'e derived as explained in Section im 
taking into account BHB star probabilities and distance uncertainties. The 
four panels illustrate the dependence of these maps on the probability cut. 
The prominent overdensities in these maps are conserved for all probability 
cuts, showing that issues with the probabihty assignment are unlikely to affect 
our analysis significantly. 

Stream, but a late-type rotating disk galaxy. This model also 
reprodu ces a bifurcation in the leading arm of the stream as 
seen by lBelokurov et alJ (l2006ah . 

3. RESULTS 

3.1. Probabilistic BHB Density Maps 

We then create maps to visualize the BHB star density, 
which account both for the finite probability that stars are 
BHB stars (as described in Section 12.1b and for the distance 
uncertainties. We account for the distance uncertainty by 
viewing each star as an ensemble of 100 sub-objects, with 
line-of-sight distances drawn from a Gaussian distribution 
with 5% scatter (see Section l2.1.3l l around the mean of the 
distance estimate for each individual BHB star. These sub- 
objects, each of which has a probability of 1 %, are then multi- 



plied by the probability they have to be a BHB star. Through- 
out this paper, we only consider stars with pbhb > 0.5 unless 
stated otherwise. A map in the Sgr plane is then created by 
dividing the plane into cells for which the probabilities of the 
included stars are summed. The 'signal' therefore depends on 
both spatial abundance of stars and on probability of each to 
be a BHB star. These maps then get convolved with a Gaus- 
sian kernel with a size of cr = 1 kpc for presentation purpose^ 
We apply this technique to create spatial maps of the Sgr de- 
bris and for plotting the velocity distribution along the orbital 
longitude of the system. In Figure |5] we illustrate the effect 
of different probability cuts on the Sgr plane. T his is of par- 
ticular interest in the context reported in Section ETI that the 
probability assignment is assumed to work less well for larger 
distances. 

Our basic map, the distribution of BHB stars in the Sgr 
stream plane, is shown in Figure |6] (top panel). The upper- 
most panel shows the full sample of stars, where the overden- 
sities are pointed out by dashed lines. Clearly visible is the 
leading arm of the stream to the right of the plot (white line). 
Less prominent, but still significant (see belowfl is the over- 
density denoted by the black dashed line. In common with 
lNewbei:g_et a l. (200 3) wh o detected part of the overdensity 
and iNewberg et al.l (l2007l) where it was also shown in BHB 
stars, we provisionally attribute this to the Sgr trailing arm. 
We find further support for this overdensity in the on-sky plot 
of a broad distance slice (60 kpc < li < 80 kpc) covering most 
of the overdensity seen in the plane. Figure |7] shows the on- 
sky view in which the plane is clearly visible as a overdense 
region. Also clearly visible in Figure |6] is the globular clus- 
ter NGC 2419 at (x,y)=(-80,-35), but its relation to the Sgr 
trailing stream is unclear. 

In the direction of the leading arm, BS contamination is 
faintly visible as an echo of the leading arm at ~ 80 kpc 
from the Sun. It is noteworthy that our selection has 
significantly reduced this con tamination compared to, e.g., 
iNiederste-Ostholt et al.l (120101) . Since we assign all stars with 
Pbhb > 0.5 as BHB stars this causes an overestimate of BS 
star distances (by 1-2 mag, or a factor of two or so in dis- 
tance, as observed). We illustrate the expected location of this 
shadow caused by stars in the leading arm region by the dotted 
black lines in Figure |6] (giving the transposition of the white 
line for stars overestimated by 1.5 and 2 mag, respectively). 
We show also the Bootes dwarf, which happens to lie in the 
Sgr plane. 

We have adopted two different methods to estimate the sig- 
nificance of the candidate trailing stream. In the first ap- 
proach, we estimate the significance in small areas of 4 kpc 
X 4° along the trailing stream. We divide the plane into areas 
of constant radial and angular extent, and count the number 
of stars in these fields. For a field / the number of stars in the 
field is A^, . The mean number of stars in a ring with constant 
heliocentric distance Nring and standard deviation of cTrmg is 
derived for each value of heliocentric distance range to ac- 
count for the increasi ng vol ume of the wedge with increasing 
distance (see Section 12.1.21 for a description of the geometry 
of the plane). We also exclude the angular range to the right 
of the area indicated in Figure |6] from the calculation of the 

^ Later we will use a polai' coordinate system, which is defined in Sec- 
tion l3.2l where a kernel of 0.5 kpc in distance and ~ 2° in the orbital angle 
coordinate is used. 

' Among other features this appears somewhat more prominent if a lower 
probability cut for BHB stars is applied (e.g., 20% or 30%, see Figure|5). 
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Fig . 6. — Density maps of stars with pbhb > 0.5 in a plane which includes the Sgr tidal tail and the Galactic Center. The maps are derived as explained in section 
I3.1l taking into account BHB star probabilities and distance uncertainties, with a minimal accepted probability of 50%. The upper panels shows the full sample 
whereas the two lower panels show the sample split into redder (right panel) and bluer (left panel) BHB stars. The Sun is located at the origin. To enhance the 
contrast in the regions of the Sgr debris the innermost 20 kpc are not shown. Clearly visible is the leading arm of the Sagittarius stream stretching approximately 
from -10 to 30 in x and —20 to —40 in y (indicated by the short white dashed line). This feature is very clear in the full sample and the red subsample whereas 
it is very patchy in the blue subsample. The two outer dotted black lines illustrate where the expected shadow of misinterpreted blue straggler stars should be 
assuming stars which are intrinsically 1.5 and 2 mag fainter, respectively. Another clearly visible feature is the cluster NGC 2419 at (x,v)=(— 80,-3 5) which is, in 
contrast to the leading arm, more prominent in the blue. It is not clear if NGC 2419 is associated with the candidate trailing arm (see Section[3]4). We also see a 
very faint indication of an overdensity in the region of —80 to 20 in x and around -50 in y, spanning the region between the cluster and the leading ami (indicated 
by the long black dashed line). If this overdensity is real it could possibly be associated with Sgr and repres ent p art of the trailing arm. The upper right panel 
illustrates the selection region used for an estimate of the significance of this overdensity discussed in Section lTll The Bootes dwarf galaxy can be clearly seen 
just above the leading arm of Sgr at (x,y)=(20,-60). 



mean and standard deviation for all distances to avoid the ob- 
vious overdensities from the leading arm in this area as well 
as the contamination at larger distances from misinterpreted 
BS stars. The significance of any deviation in the number 
of stars of each field within this sample of equidistant fields 
in units of the standard deviation a for region / is given by 

ai = (Ni-Nring)/(Tri„g. 

We take the region around the suggested position of the 
trailing arm as indicated in the upper right panel of Figure 
|6]by the dash-dotted lines and compare the average deviation 
of these fields with a comparison sample in the same plane 
but outside the trailing arm area. Note that this area does not 
include NGC 2419 to get a clean estimate of the significance 
of the proposed trailing arm. These fields are chosen in a way 
that the number of fields per distance interval of on- and off- 
stream fields is the same. The 57 on-stream fields show an 
average deviation of +QAa per field, indicating a weak over- 
density, whereas the 57 off-stream fields show with an average 
deviation of -0.6(t per field the corresponding underdensity. 

To get an idea of the significance of the whole extent of the 
structure we adopt a larger area, as shown in Figure (up- 
per right panel). Within this region consisting of 200 fields 
we randomly select a number of fields, equal to the number 
of stream fields we used earlier, and determine the average 
number of stars in this selection. Applied many times this 
bootstrapping method gives an estimate of the mean value and 



standard deviation we can expect in a randomly selected struc- 
ture of this size. We find a mean value of 23.5 stars per field 
in the large box with a standard deviation of 1.3 stars. The 
average number of stars in the selected structure fields is 28 .5 
per field which corresponds to an deviation of 3.8 a from the 
mean value. The candidate stream fields are compared with 
all fields - including stream fields - potentially underestimat- 
ing the significance. 

In the two lower panels of Figure |6] we show the maps that 
result after splitting the BHB sample in g-r color at g- r = 
-0.18, such that the number of stars with g < 18.5 is about 
equal in the red and blue subsamples. The main motivation 
to do so is to probe the variations of the stellar population in 
the Sgr stream. In Figure |6] we show the red subsample in 
the lower right panel and the blue subsample in the lower left 
panel. We find that Sagittarius (especially the leading arm) 
is much more prominent in the red stars (see Figure while 
other parts, such as NGC 2419, are dominated by blue stars. 

In summary, we find Sgr's leading arm to be a prominent 
feature in BHB stars, even more so when the BHB star sam- 
ple gets restricted to stars which are on the red part of the blue 
horizontal branch in g- r color Furthermore we observe a 
faint overdensity stretching out over most of the plane cov- 
ered by SDSS, connecting the leading arm with the globular 
cluster NGC 2419. This ov erdensity was also described by 
iNewberg et aH (120031 l2007h as a pan of the ti-ailing arm of 
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Fig. 7. — On-sky view of a distance shell (60 kpc < d < 80 kpc) which 
holds most part of the faint overdensity we see in the plane which could 
represent part of the trailing arm of Sgr. The borders of the plane are indicated 
by the dotted lines, for better visibility they are indicated by thick solid hnes 
outside the survey field. The stars that we attribute to the Sgr trailing arm 
with 1 80° <l < 270° (and what we argue are blue stragglers in the leading 
arm at 270° < / < 360° ) are cleai'ly confined to the Sgr plane. 

Sgr. 

3.2. Thickness of the Leading Arm, and Spatial Selection of 
Sgr BHB Star Candidates 

In this section we measure the line-of-sight thickness of the 
Sgr leading arm and use this measurement to select a sam- 
ple of highly likely Sgr member stars. In this subsection, we 
will present a selection based on the spatial distribution only 
which will be used for the analysis in the following subsec- 
tion. Later we will restrict the selection of a 'clean sample' 
to the radial velocity subsample for which we apply a similar 
selection technique. In what follows, we restrict our attention 
to the Sgr leading arm; the trailing arm (and candidate trailing 
arm debris) is in the wrong hemisphere and/or too distant to 
have SDSS radial velocity information. 

We adopt the heliocentric polar coordinate system defined 
by Maiewsk i et al. (2003) which was also used by Law et al. 
(|2005) and iLaw & Majewskil (1201 Oah . In this system, the an- 
gle is defined as Aq = 0° passing through the main body of 
Sgr and increasing along the direction of the trailing tail of 
Sgr The definition of the coordinate system is illustrated in 
Figure |8] where the prolate version of the models is shown to- 
gether with our BHB star data in the large panel. The inset 
panel shows the data alone. 

We first measure the width of the leading arm using the full 
sample of stars in the Sgr plane (not the Gaussian-distributed 
sub-objects), as is appropriate for measuring line-of-sight dis- 
tance scatter; the distribution of the stars is shown in Fig- 
ure |9] We divide the angle-distance-plane into angular slices 
along an orbital angle range of 250° < A0 < 300° and fit 
the distance distribution of the stellar density with a function 
consisting of three components: an exponential function and 
constant component to fit the background distribution of halo 
stars, and a Gaussian for the Sagittarius stream. The fit can be 
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Fig. 8. — Compaiison of the BHB star density (> 50% sample) with the 
Law et al. (2005) prolate potential model of the Sagittarius stellar stream and 
in the small panel data alone for clarity. The colors showing the model parti- 
cles show debris lost from the progenitor during different epochs. The yellow 
points show debris stripped since the last apogalacticon, while green, blue, 
and purple show debris which became unbound two, three, and four orbits 
ago, respectively. The Sun is located in the center of the coordinate system. 
The dash-dotted line gives the position of the Galactic Plane, with the orien- 
tation is chosen such that it falls on the x-axis. The solid black line indicates 
the direction to the Sgr dwarf galaxy which also defined the zero direction of 
the longitudinal coordinate system, with the angle Aq increasing clockwise. 
To enhance the contrast in the regions of the Sgr debris, the innermost 20 kpc 
are not shown. Except for the 'leading arm', any data-model correspondence 
is not obvious. 

described by the expressioifl S((io | A0) = P()+P\ * exp(-^)-l- 

Pi * exp(-0.5 * i^^f). The best fit was determined using 
a chi-square algorithm. As the Sagittarius leading arm is sig- 
nificantly more prominent in the red subsample of BHB stars 
we also apply the fit to the red part alone (see Figure [TOt. For 
comparison the histogram of the corresponding distribution in 
the models (with a prolate potential) is shown by the dashed- 
dotted line. The histogram is scaled down by a factor of seven 
to approximately match the number of stars in the data. The 
models show a bifurcation of the leading arm in distance be- 
tween the debris lost in different orbits. We cannot see this 
in our data, the relative separation and size of the peaks are 
roughly of the same size as the fluctuations we see in the data 
in a typical angle slice. The results of the Gaussian fit are 
shown in Table|2]and Figure|9j crosses denote the mean value 
and the 'error' bars show the standard deviation a around that 
mean. 

We use these results as a first step in isolating a clean sam- 
ple of BHB stars. A second-order polynomial is fit to the 
mean values, shown in Figure |9] We use this line, shifted 
by ±2 times the mean standard deviation as borders w ithin 
which we select leading arm member stars. In Section 13.51 
we will refine this selection by taking into account an addi- 
tional selection in velocity space. In the following we will use 
the spatially selected sample defined here since the kinematic 
selection also very strongly limits the sample size to stars 

■* Due to the proximity of Bootes, we added a second Gaussian to this ex- 
pression to isolate the profile of Sagittarius in the relevant bins from Bootes. 
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Fig. 9. — Heliocentric distance as a function of orbital longitude (see also the left column of Figure [13) ■ In the left panel, the prolate model is shown as gray 
points. In the right panel, the data are shown with the small green points representing the photometric sample (limited to BHB probabihties > 50 %) a nd the lai'ger 
orange points representing the BHB sub-sample with radial velocities. The result of the Gaussian fit to the red BHB stars in the sample (Figure [To) is shown as 
asterisks (mean distance) and 'error bars' (width a). We fit the mean values with a second order polynomial function (central line) and take ±2(7 as our distance 
cuts (outer two lines). The dashed vertical lines mark the ran ge in Aq over which the fit to the angle slices was performed. The orange dots falling within these 
outlines we denote as the kinematic BHB selection for Figure[T4l 



TABLE 2 
Width of the Leading Arm 



(Aq) 


{d) prob>0.5 


T pmb>0.5 


^ intr,pi'oh>0.5 


{d)red 


Cred 




(degree) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


252 


36.5 


2.8 


2.1 


38.6 


1.2 




258 


35.6 


1.3 




34.7 


2.0 




264 


42.9 


4.4 


3.9 


43.5 


2.3 




270 


46.7 


6.2 


5.7 


44.3 


5.0 


4.2 


276 


46.5 


3.6 


2.8 


47.6 


4.1 


3.0 


282 


48.5 


5.0 


4.4 


48.5 


4.8 


3.7 


288 


50.2 


4.8 


4.1 


50.4 


3.8 


1.7 


294 


52.7 


6.5 


6.0 


53.1 


5.6 


4.4 



Note. — Results for the Gaussian fit to the leading arm. Mean (d) and standard deviation a values are given both for the 50% probability sample and for the 



red subsample. We present also an estimate for the intrinsic width ai„tr = sqrt{cr^i^^- {0.05(d))- -cr^ 



bootstrap 



which have radial velocity data available. The spatially se- 
lected sample will be limited to BHB star probabilities greater 
than 50%, whereas no probability cut is applied for the radial 
velocity sample since these stars are spectroscopically classi- 
fied. 



3.3. 



Bifurcation of the Leading Arm Perpendicular to the 
Plane 



Following lYannv et"an (120091) . we also look into the 
bifurcation of the lead ing arm as it was discovered by 
iBelokurov et al.l(l2006 aV When looking at thin distance slices 
of the SDSS using a population with a high abundance like 
MSTO stars one easily sees that the Sgr stream splits up into 
two parts. Given the relatively sparse distribution of BHB 
stars in the Sgr stream we use MSTO stars to define a se- 
lection for the two parts (see Figure [TH . In our BHB star 
sample its elf we do not se e any indication for a bifurcation. 
Following Bell et alJ (I2008 ). we select MSTO stars in a color 
range of 0.2 < g-r < 0.4 and a distance modulus range of 
16.5 < m-M < 17.5 assuming an absolute magnitude of 
Mr = 4.5. This corresponds to a distance range of 20 - 32 kpc. 
In Figure [TT] we show the distribution of these MSTO stars. 
In the upper panels, we also show the distribution of the Sgr 
BHB stars as spatially selected from Figure |9] Note that the 



clearly identifiable part of the leading arm in BHB stars is not 
in the region on the sky where the bifurcation is most appar- 
ent. In the lower right panel, we illustrate the low density of 
BHB stars in the relevant distance slice (same distance modu- 
lus selection as for MSTO stars) with A0 < 250°, which pre- 
vents us from investigating this part of the leading arm in BHB 
stars. Consequently, in the following we only study the lead- 
ing arm for Aq > 250°. In Figure [12] we present measure- 
ments of the mean and width of the two branches in thin angle 
slices. Note that in contrast to Figure|9]this measurement was 
made on the pre-selected sample and not fitted to the data in 
the same fashion as illustrated in Figure[TO] The two branches 
show similar distances with a 1-2 kpc variation in the mean 
distance values (see also Table [3] for a listing of the results). 
Several studies showed a systematic separation in the distance 
of the two branches, such that the high galactic latitude part 
of the stream is closer for most of the leading arm as seen in 
the SDSS. Y anny et all (12009) report this offset in the distance 
distribution of BHB stars along the leading a rm by visual im- 
pression. The same tre nd was also seen by IB elokuro v et al.l 
(2006a) (results Hsted in lNiederste-Ostholt et al. 2010), show- 
ing an offset of 2-3 kpc. An offset was also given by the 
iPenarrubia et al.l(l2010b ) models for which we show the mean 
distances of the two branches in Figure [12] separated and 
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distance (kpc) 

Fig. 10. — Distribution of heliocentric distances for 6° wide angle slice, 
285° < Aq < 291°, of Sgr's leading arm. The small panel on the upper right 
illustrates the angle slice used for this histogram. The full sample is denoted 
by the empty solid line histogram, while the gray histogram represents the 
sample with a >50% probabihty cut. A Gaussian has been fit to represent 
the stream stars in this subsample, while the fore-/back-ground has been fit 
by a power-law; the resulting distance and line-of-sight thickness are noted 
in the top left. Shown in the histogram filled with red vertical lines is the 
subsample of red BHB stars and shown in the one filled with blue horizontal 
lines is the subsample of blue BHB stars. The leading ann shows a slightly 
larger number of stars from the red subsample for all angles. The dash-dotted 
line shows the histogram for the prolate Law et al. (2005) model in the same 
angle slice. It has been scaled down to match approximately the number of 
objects in the data. 

measured in the same way as our data. Although we do not 
see a clear separation in distances in our data, the mean dis- 
tances of the two branches and their relation to each other are 
sensitive to small change s in the separation c ut between the 
two branches. Recently, Correnti et al.l (12010) measured the 
distances of these two branches in Red Clump stars finding 
also only a small offset between them which is of a similar 
order as found in this study (see, e.g., their Figure 13). 

3.4. Kinematics and Comparison to Models 

To improve our understanding of the origin of overdense 
regions of the Sgr plane and the likelihood of those over- 
densities being associated with the Sgr system, we compare 
our data with models of t he Sg r debris by La w et al. (2005), 
iLaw & Maiewskil ( l2010ah . and lPenarrubia et al„ C2010b ). To 
complement our dataset with kinematic information, we use a 
sample of radial velocities determined from SDSS DR7 using 

TABLE 3 

Mean values with uncertainties and width for the two 

BRANCHES OF THE BIFURCATED STREAM AS PRESENTED IN FIGURe[T21 



(Aq) 


{<^)high 




{d)iow 




(degree) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


252 


374 ± 1.4 


3.9 


36.8 ± 1.2 


3.1 


258 


364 ± 1.9 


4.0 


37.2 ± 1.5 


34 


264 


39.9 ± 0.8 


3.9 


4L4± 1.1 


4.6 


270 


42.5 ± 0.8 


3.6 


42.9 ± 1.1 


4.7 


276 


45.2 ± 0.7 


3.9 


47.6 ± 0.9 


3.5 


282 


48.4 ± 0.9 


4.3 


484 ± 14 


3.4 


288 


51.0 ±0.6 


3.7 


50.3 ± L4 


3.8 



the method of IXue et alj (l2008h and lXue et all (1201 Oh (for an 
illustration of the spatial distribution of these stars with kine- 
matic information in the Sgr plane, see Figures[T3]and|9). The 
models are based on the 2MASS data set of M giants, see Fig- 
ure|4]for their distribution. Figure [T3]shows probability maps 
made in the same fashion as described above for both distance 
and velocity as a function of the orbital angle. In the thre e 
lower panels, we overplot the models by iLaw et akl (120051) . 
From top to bottom, they show the tidal debris in a prolate, 
spherical, and oblate Galactic halo potential. In the third and 
fourth rows, two models using a triaxial halo potentia l are 
shown; in the fourth row t he model from [L aw & Ma iewskil 
( I2010al) and the model by iPefiarmbia et alT (,2010b) "in the 
third row which, in contrast to the other models, use a disk 
galaxy as a progenitor of the Sgr stellar stream. In th e sec- 
ond ro w, a radial velocity sample of M giants (Majewsk i et al.l 
120041) is shown alongside with our BHB star data set. They 
cover mostly parts of the stream not covered by the SDSS. 

In the following, we compare the location of the predicted 
debris in the different models with our observations of the 
distribution of BHB stars. This comparison is merely meant 
to illustrate tentative agreements and disagreements between 
the models and the data with the goal of identifying features 
and regions of interest for further investigation, and not give 
a conclusive answer for a best model. 

The prolate and triaxial models, which are shown in the 
third to fifth rows of Figure [13] clearly show the best con- 
sistency with the leading arm which is the most prominent 
part in the SDSS BHB sample (at an orbital angle of 230° < 
Aq < 300° and heliocentric distance betwe en 20 and 60 kpc) . 
On the other hand, the trailing arm from the iLaw et ail (12005b 
spherical and oblate models stretches out to larger distances 
than the prolate model, qualitatively (but not quantitatively) 
matching better the candidate Sgr debris shown in Figure lafl. 
The recent models by Peiiarmbia et al. (2010b) show a trail- 
ing arm which stretches out to much larger distances than in 
the other models. Still, w e do not see a good match with 
the observed overdensity. ICorrenti et al.l (12010b report de- 
tection of a trailing arm segment in Red Clump stars which 
appears to be consistent with the prolate models around the 
crossing region of the leading and trailing arm in the range of 
220° < Aq < 290°. This feature is observed at much smaller 
distances than what is suggested here. We do not focus on this 
distance range here, as in at least our investigation we find a 
high degree of contamination from a nd/or cross-talk with t he 
Virgo overdensity. If the detection of ICorrenti et al.l (120101) is 
interpreted correctly as part of the trailing arm the overdensity 
seen here could possibly belong to a different trailing wrap. 

Turning to the possible association of NGC 2419 with 
the candidate trailing arm debris, we note that the he- 
liocentric radial velocity of NGC 2419 was measured 
by iPeterson et al.l ([1986) to be -20kms"' which corre- 
sponds to a galactic standard of rest velocity of -14km s"' 
dNewberg et al.l l2003b . This corresponds well with the hy- 
pothesis that the c luster is a part of the trailing stream near 
its apogalacticon (iNewberg et al.l 120031) . In addition, the 
properties of NGC 2419 are unusual in its own right (e.g., 
iDalessandro et al.l 120081) : it is very luminous with My ^ 

^ Part of the issue in reproducing such debris may be related to the dis- 
tance offsets between the BHB stars and M giant tracers of the Sgr tail. The 
models were built to reproduce the smaller distances characteristic of the M 
giant tracers; we speculate that models reproducing better the leading arm in 
BHB stars would more easily yield a trailing arm consistent with the distant 
candidate Sgr debris. 
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Fig. 1 1 . — Map of the MSTO stars in a thick distance modulus slice, 16.5 < m—M < 17.5, corresponding to a distance range of 20 — 32 kpc (left side and top 
right comer). The dotted black lines indicate the borders of the plane as this projection represents a side-view of this plane. The orbital angle Aq is given on the 
right of the plane segment. In the top panels, the BHB stars which were spatially selected as member stars of the leading arm are overplotted in yellow and blue. 
For better visibility, the area of interest (red box) is enlarged in the top right panel. The two colors denote the selection for the two parts of the arm following the 
appearance in MSTO stars. In the bottom right panel the BHB star distribution in the same distance shell is shown to illustrate why we can not select stars in the 
closer part of the stream where we see the strongest bifurcation in MSTO stars, as the BHB stars get much less abundant in the corresponding angle range (lower 
half of the plot). 



-9.5 and has a large half-light radius r/, ^ 25 pc jBellazzinil 
120071) . placing it in a region of radius-luminosity param- 
eter space populated also by wCen and M54, that have 
both bee n argued to be the stripped cores of dwarf galax- 
ies (e.g.,|Sa raiedini & Layden 1995; Hilker & Richtler 2000^, 
Romano et al. 2007; Bellazzini et al. 2008; Georgiev et al. 
2009I) . Yet, the situation with NGC 2419 in particular is not 
clear cut. There is no evidence of m ultiple stellar popula- 
tions in NGC 2419 jCohen et al.ll2010l) . i n apparent contrast 
with the properties of, e.g., wCen (e.g., Ripepi et al. '2007t 

tndquist & Hess 2008). Furthermore, Casetti-Dinescu et al] 
)09) have calculated a preliminary orbit for the Virgo stellar 
overdensity, finding that it is very eccentric, and they suggest 
that NGC 2419 may in fact be associated with the Virgo stel- 
lar overdensity rather than Sgr Furthermore, we do not see 
a clear velocity signature of trailing debris in the SDSS ve- 
locities (although it is unclear if a signature is expected in the 
sparsely-sampled SPSS BHB veloci t y data set). Finally, the 
updated models of iLaw & Maiewskil ( l2010al) in a triaxial po- 
tential show an increased inconsistenc y with NGC 2419 as 
described by ILaw & Maiewskil (I2010bl) . 

Although the full velocity sample as we show it in this plot 
does not show a very clear signal for the prominent leading 
arm, it is still obvious that the main overdensities (240° < 



A0 < 300°) agree best with the models for the prolate and 
triaxial versions. This will become clearer when we restrict 
the velocity sample to stars within the region of the leading 
arm in distance space in the next section. 

A serious inconsistency with the models can be seen in the 
region where the trailing arm is predicted to stretch into the 
region covered by the SDSS (around (-60,0) and upwards in 
Figure [8). In the data we do not see a signal which would 
come anywhere near the intensity which is predicted by the 
models for this part of the arm. The absence of such a coun- 
terpart indicates a serious problem with the models. This can 
not be explained through differences in the stellar populations 
in the debris; the models predict this part to consist of stars 
that got unbound in the same orbits as the debris in the part 
of the leading arm that we can observe in the SDSS. We spec- 
ulate that this discrepancy may be alleviated in models tuned 
to reproduce better the distances of the leading arm as traced 
by the BHB stars. 

In the following section, we attempt to measure the veloc- 
ity spread of the Sagittarius stellar stream. We continue our 
attempt to isolate a 'clean' sample of stars most likely belong- 
ing to the Sagittarius stellar stream. We use both positions and 
kinematics to achieve a high reliability of our selection. How- 
ever, the size and distribution of the radial velocity sample 
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Fig. 12. — Mean distance and distance error corrected line-of-sight thick- 
ness of the two branches as determined from the spatially selected sample 
of Sgr BHB stars. Blue represents the branch at higher ga lacti c latitude and 
yellow the one at lower galactic latitude (see also Figure II U . For presen- 
tation purposes, the points are offset by half a degree to the left and right, 
respectively. The uncertainty of the mean value determination is determined 
by bootstrapping and given by the solid errorbars. The dashed 'errorbars' 
represent the width of each substream, also corrected for the sample selec- 
tion effects via bootstrapping. For better visibility, they are offset by another 
half a deg ree. The lines show the results of the same measurement on the 
[Pefiarrubia et al. (2010b) models. We do not see a clear trend in the distance 
offset of the two branches as in the models and the observed offset is much 
smaller than predicted by the models and seen in other analyses with different 
stellar populations (e.g., Niederste-Ostholt et al. 2010). 

limit this selection strongly. 

3.5. Selection of a 'Clean Sample' of BHB Star Candidates 

In this section, we continue our effort to select a 'clean sam- 
ple' of Sgr BHB stars. In Section [J!2l we already made a spa- 
tial selection of the leading arm stars. In the following, we 
will restrict this selection to the radial velocity subsample to 
achieve a sample which follows the leading arm in both dis- 
tance and velocity space. This selection of a 'clean sample' 
of Sgr BHB stars is based purely on the data, but agrees qual- 
itatively in both distance and velocity space with the Law et 
al. models. 

Figure [14] shows the radial velocity full sample (orange), 
with those lying in the distance selection in red. As can be 
seen in comparison with the models (in gray), the selected 
stars are mostly concentrated in an area quite consistent with 

TABLE 4 

Radial Velocities of the Leading Arm 



(A0> 
(degree) 


(RVgal) 
(km/s) 


O'obs 

(km/s) 


(km/s) 


254 


-76.2 


40.7 


35.3 


264 


-28.4 


100.7 


98.3 


274 


-27.6 


92.8 


88.1 


284 


-5.7 


37.3 


35.8 


294 


24.6 


43.3 


40.2 



Note. — Results for the Gaussian fit to the radial velocity distribution of 
the leading arm sample. The intrinsic velocity spread is corrected for the ra- 
dial velocity error and the sample composition in the bin (via bootstrapping). 



the general trend of the model. We isolate candidate leading 
arm stars by taking angle slices in which we fit a Gaussian 
to the distribution of the distance-selected stars (red). Again 
we fit the mean values with a second-order polynomial func- 
tion (see Table H] for results). In Figure [15] we show the 
trends in distance and velocity, and their intrinsic dispersions, 
compared with the models with prolate and triaxial potentials. 
Stars lying within the distance selection, and within one mean 
standard deviation in either direction of the velocity fit are 
included in our 'clean' sample (see Table [5]for a full listing 
of the objects). When using this sample for further analysis 
one has to keep in mind that it is strongly restricted by the 
uneven coverage and magnitude distribution of the radial ve- 
locity sample. 

4. DISCUSSION AND CONCLUSIONS 

A number of previous works have explored the properties 
of the Sgr tidal debris using different stellar tracers, such as 
main-sequence tumoff stars, M giants, subgiant stars, or BHB 
stars. In this paper, we have presented a discussion of the 
structure and properties of the Sagittarius stellar stream us- 
ing candidate BHB stars selected from the SDSS coverage of 
the North Galactic Cap. BHB stars are in many senses an 
excellent tracer of tidal structure: they are luminous and can 
be traced to > lOOkpc distances from the Sun with current 
surveys; they are good standard candles with ~ 5% accurate 
distances; and although they are rather sparse compared to 
other stellar populations, they are still quite numerous in the 
Sgr tidal stream (with exception of the closest part of the lead- 
ing arm as shown in Figure [TTJ. Recently, there h ave been a 
number of Sgr stellar stream studies (e.g. , Yannv et al.ll2009l ; 
iNiederste-Ostholt et al.ll2010l;ICorrenti et al..,2 010) also based 
on various stellar populations, some including BHB stars in 
SDSS. 

In contrast to those studies, we entirely focus on BHB stars 
for the analysis, using other data only to relate our distance 
scale to other distance indicators. We attempt to use as pure 
a sample of BHB star candidates as is possible for our anal- 
ysis, minimizing to the greatest extent possible the high lev- 
els of contamination seen in earlier studies. For charting out 
the global structure of the Sgr tidal stream we make use of a 
method that selects BHB stars from SDSS imaging data us- 
ing a spectroscopic training set to isolate areas of ugr color 
space that give a sample that should consist of ^ 70% BHB 
stars. This method does not just make binary acceptance or 
rejection decisions based on the position in color-space, but 
assign probabilities to the stars based on their position in color 
space. In our analysis we mostly reject stars with probabilities 
< 50% and make use of the probability information for the 
remaining sample by weighting the individual stars by their 
probabilities. 

We evaluate the precision of our distance determination 
through comparison with distance measurements of known 
clusters and dwarf spheroidals. We see an offset in the mean 
values of ^ 4% for the literature values, and a distance vari- 
ance of '^5%. Comparison to the M giants in the Sgr orbital 
plane implies that the M giant distances should be revised up- 
ward by 8% or 12% when compared with our BHB scale or 
the cluster distance scale in the literature, respectively. The 
offset to the previously adopted M giant distance scale is ex- 
pected to propagate throu gh to the models built to match the 
M giant observations (e.g.. lLaw et al.ll2005tlLaw & Maiewskil 
[2010ai) . When studying the kinematics of the Sgr tidal tail, we 
focus on a sample of stars with SDSS DR7 spectroscopy clas- 
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Fig. 13. — Distance and velocity of the Sgr stream for data and models as functions of orbital longitude. The upper row shows the data alone for stars within 
the Sgr orbital plane, analogous to Figure |6] The right column shows the distiibution of measured heliocentric radial velocity in the Galactic standard of rest for 
the small subsample of BHB stars with SPSS spectroscopy. T hese data are repeated in the background of each row. In the second row from the top we show a 
sample of M giants with radial velocities froml Maiewski et aTl 12004 ) in red (right, and the spatial distribution on the left), while we show the spatial distribution 
of our BHB star radial velocity sample in orange on the left. In the same row, the velocity of NGC 2419 is represented by the dark red star In the lower five rows, 
we show different models of the Sgr debris. In the third and fourth row s , the t riaxial halo models by Law & Maiewski (2010a) and Pei iarrubia et a l. (2010b) are 
shown. The fifth through seventh panels show the models b v l Law et alj j2003h using different halo potentials: prolate, spherical, and oblate. The different colors 
in the models show debris from different orbits (see Figure |8). Comparing the BHB star maps with the models we see the best match for the leading arm region 
in the triaxial and prolate models for the spatial distribution. Although the radial velocity map does not resemble very well what is expected by the models, the 
main overdensities are also best covered by the triaxial and prolate models. Apart from the leading arm there is no good match between any of the models and 
the data. The overdensities in BHB stars seen at distances between ~ 50 and 90 kpc is not reproduced quantitatively by the models, nor do we see overdensities 
in the regions the models predict for the trailing arm. 



sified as BHB sta rs using the method of lSirko et alJ (120041) and 
iXue etal.l(l2008h : this sample should be > 90% BHB stars. 

With these samples, we focus on four Sgr stream issues that 
are not well-explored in the literature: a possible extension of 
the trailing Sgr debris stream, the line-of-sight thickness of 
the leading tail, the bifurcation of the leading arm and the he- 



liocentric distances of the two branches in BHB stars and the 
isolation of a small sample of high-probability Sgr member 
stars. 

Using the photometric sample with a > 50% BHB probabil- 
ity limitation, we identify a possible extension to the trailing 
tail of the Sgr debris stream to 60-80 kpc. The densest part of 
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Fig. 14. — Heliocentric radial velocity in the Galactic standard of rest as a function of orbital longitude for the 'Law et al.l (120051) prolate model and for BHB 
stars in the leading arm. The radial velocity sample of BHB stars is shown by the orange and red dots (right panel); big red symbols denote stars falling into the 
distance selection box of Figure |9] and small orange denotes everything falling outside this selection region. The majority of red dots follow roughly the trend 
indicated by the models (left panel). We fit the velocity distribution of red symbols with a Gaussian, using somewhat wider angle slices than before to account for 
the sparser sample. The mean and standard deviation of these fits are given by the asterisks and the bars (see also Table|4]. We apply a second order polynomial 
fit to the mean values and shift this by the mean standard deviation in both directions to identify a selection ai'ea for the 'clean' sample. 




Fig. 15. — Distance and velocity dispersions of data from the spatial and kinematic selecti ons for the leading arm and models. The veloc ities are again given 
in the Galactic standard of rest. We only show the models for a prolate and triaxial potential l IMaiewski et al.l2004HLaw & Maiewskill2010al) since these showed 
the best consistency with the properties of the leading arm. The data points represent the mean and width measurements (bars) as presented in Figures ID and 
1141 showing a fit to the total sample on the left and the spatially selected sub-sample on the right. The data are corrected for distance uncertainties and for 
uncertainties resulting from the sample selection in the bins by bootstrapping. For some bins these uncertainties are too large to resolve the intrinsic width of the 
stream. 



this feature, which coincides spatially with the globular clus- 
ter NGC 2419, was pr eviously argued to be associated with 
the Sgr traihng arm bv lNewberg et al.l (120031) . Our BHB star 
maps confirm a weak overdensity which may be the extension 
of this arm back towards the Milky Way, which was also seen 
by Newberg et al. (2003, 2007). We estimate the significance 
of this feature to be around 3.8 a as compared to random se- 
lections of the same area within a region spanning the angular 
range and distance of the proposed trailing arm. A concentra- 
tion in this region, which is claimed to be associated with the 
trailing arm of Sgr was also found by Sharma et al. (2010) in 
2MASS M giants through a group finding technique. Such a 
feature is expected qualitatively by models of Sgr disruption. 
Quantitatively, all models predict that this 'returning' segment 
of the trailing arm should be closer to the Sun along these lines 
of sight. Yet, BHB stars are not observed at these predicted 
distances, and this tension would be resolved if one instead in- 



terpreted this distant overdensity as this predicted part of the 
trailing Sgr arm and acknowledges a discrepancy between the 
positions predicted by the models and the obs erved location. 
In this context it is worth noting that recently ICorrenti et al.l 
(1201 Oh reported an overdensity in Red Clump stars consistent 
with the predicted trailing arm location in the prolate models 
in a range of 220° < A0 < 290°. Owing to confusion be- 
tween Sgr and Virgo overdensity debris at the distance ranges 
probed by ICorrenti et al.l (120101) . we were unable to confirm 
or refute this feature. If their feature is indeed correctly in- 
terpreted as trailing arm debris, we would suggest the feature 
identified here may be another, more distant wrap of trailing 
arm debris from an earlier close passage of Sgr 

We use the > 50% probability sample to characterize the 
leading arm of the Sgr stellar stream more closely by measur- 
ing the line-of-sight thickness and selecting a high-probability 
sample of member stars. We find a mean thickness of ~ 3 



The Structure of Sgr in BHB Stars 



TABLE 5 

Sample selected for a high probability of association with Sagittarius. 



R.A. 


Decl. 


/ 


h 


g 


u-g g-r pbhb 


HRV HRVerr 


RVgal Fe/H Fe/H error 


(deg) 


(deg) 


(deg) 


(deg) 


(mag) 


(mag) (mag) 


(kms^l) (kms^l) 


(km s"l) 



149.31784 
149.32185 
151.66303 
151.68573 
152.46681 
152.73950 
162.91279 
163.11225 
165.03563 
165.26730 
171.35795 
177.49904 
181.41167 
181.54416 
181.57671 
185.50090 
185.50929 
185.13369 
173.32034 
191.92214 
192.82045 
196.78478 
249.77752 
236.23186 
247.67761 
195.51961 
198.27333 
203.08792 
203.46564 
235.28992 
237.33833 
183.15324 
184.33572 
186.20173 
230.36425 
231.31942 
247.71507 
177.38054 
213.82899 
213.89940 
214.76100 
12.739674 
21.239977 
11.571439 
37.994698 
62.479321 
40.554837 
42.311590 
181.39053 
199.43635 
170.96747 
192.16892 
174.28068 
181.68676 
184.50026 
185.92025 
188.14379 
187.30791 
188.69635 
179.19443 
179.27744 
182.51252 
187.04228 
188.39442 
190.63810 
192.20784 
172.36979 
126.14404 
123.44799 
140.29301 



0.64158908 
0.82575434 
0.71303259 
0.83859047 
0.72481130 
0.80837915 
0.63505894 
0.80012300 
0.82406409 
0.67366782 
0.77705579 
0.69659886 
0.68832934 
0.66210870 
0.68079996 
0.76998992 
0.83180709 
0.79081299 
1.2211629 
1.1856847 
1.1740127 
1.0952469 
-0.12263493 
0.32725518 
0.34199587 
-0.73406591 
-0.68782645 
-0.74637791 
-0.74396165 
-0.66049236 
-0.64090826 
-0.35575570 
-0.30121489 
-0.34449501 
-0.33667675 
-0.34046227 
-0.35833233 
0.15937633 
0.048201690 
0.18624951 
0.13294046 
15.849297 
14.300256 
15.267388 
-9.4446953 
-6.3200403 
-8.7054956 
-8.4109077 
-3.3742984 
-2.8485374 
-2.5415740 
-2.2424141 
-3.2210852 
-3.2304553 
-2.7444692 
-2.8228684 
-2.8617039 
-2.7770480 
-2.7577484 
-2.3958845 
-2.4873776 
-2.2918592 
-2.3998936 
-2.3430983 
-2.4297249 
-2.4226168 
-1.4338628 
46.955835 
46.625513 
58.048908 



337.25807 
340.25544 
332.93978 
336.16948 
354.09098 
355.95736 
352.51890 
355.61411 
291.64879 
294.87797 
334.88970 
355.65421 
292.43919 
314.60980 
316.52701 
332.59194 
333.63472 
342.06521 
359.16254 
353.74592 
342.24188 
353.51720 
343.67858 
349.53368 
297.76277 
258.01091 
254.48333 
261.17467 
294.52133 
338.03357 
336.47791 
341.09162 
359.67778 
338.42710 
312.91600 
350.71450 
303.66505 
281.06974 
284.98406 
300.73988 
319.00435 
324.89192 
293.68604 
322.21782 
322.70253 
256.52080 
260.66485 
260.92325 
279.90974 
298.10400 
302.82768 
273.75164 
335.45465 
301.78123 
261.08378 
245.46876 
245.13008 
249.11382 
259.84128 
351.97783 
353.19843 
326.21682 
335.66563 
306.23529 
315.32404 
358.57367 
272.87946 
339.48188 
334.21365 
351.40623 



59.124824 
57.966839 
60.252136 
59.330910 
51.794464 
50.489586 
50.663479 
48.732353 
62.180853 
62.247052 
58.104219 
49.925343 
63.422670 
63.437217 
63.065209 
60.542340 
60.321873 
57.760352 
50.629827 
52.452357 
58.467647 
55.082748 
60.877749 
58.140268 
68.420754 
71.848421 
72.188679 
74.794533 
78.208821 
68.100269 
68.325007 
68.226778 
59.904530 
69.530976 
77.668905 
70.209927 
76.708831 
77.423385 
76.767400 
77.134295 
77.535698 
85.302086 
81.423698 
83.117026 
79.105115 
75.094461 
77.113803 
78.271266 
67.307502 
74.487378 
75.868349 
68.019640 
68.964352 
73.471968 
79.380844 
77.417677 
77.769645 
79.733714 
80.889018 
50.812086 
50.061905 
60.272876 
58.430511 
62.775039 
62.568188 
48.349971 
68.009004 
55.003659 
57.176227 
50.936195 



18.92 
19.20 
18.70 
19.52 
19.15 
19.11 
19.12 
19.11 
18.84 
18.87 
18.61 
18.97 
18.86 
19.01 
18.76 
18.72 
18.94 
19.08 
19.09 
19.20 
19.08 
19.13 
18.87 
19.17 
18.43 
18.38 
18.36 
18.31 
18.70 
18.99 
19.21 
18.43 
18.98 
18.77 
18.65 
18.50 
18.62 
18.41 
18.51 
18.58 
18.52 
18.74 
18.51 
18.43 
18.51 
18.11 
17.92 
18.15 
18.24 
18.30 
18.31 
18.32 
19.03 
18.58 
18.25 
17.96 
18.28 
18.85 
18.16 
19.50 
19.38 
18.84 
19.04 
18.81 
19.04 
18.96 
18.48 
18.87 
18.67 
19.12 



1.25 


-0.09 


0.68 


-5.1 


11.0 


-38.9 


-1.7 


0.2 


1.19 


-0.15 


0.72 


17.6 


19.0 


-11.7 


-10.0 


-10.0 


1.18 


-0.06 


0.44 


65.9 


11.0 


25.7 


-1.4 


0.1 


1.05 


-0.30 


0.75 


48.5 


24.0 


12.9 
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1.21 


-0.12 


0.70 


32.8 


8.0 
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0.1 
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kpc, after accounting for distance uncertainties in the BHB 
starfl comparable to the projected width of Sgr on the sky. 
These m easurements are in a similar range as those given by 
iCorrent i et al. (2010) in Red Clump stars (their Table 2), when 
the assumed overestimation by a factor of > 2 which is intro- 
duced by their measurement method is taken into account. In- 
spired by the clear appearance of the leading arm in position 
(and velocity), we use this measurement of the line-of-sight 
thickness to select a sample of highly likely stream stars from 
the spectroscopic SDSS sample. We choose stars within 2a of 
the stream in line-of-sight distance. This subsample of stars 
shows a clear overdensity in velocity space, which matches 
model predictions reasonab ly well. This strengthens the re- 
sults of lYannv et al.l (120091) who showed that BHB star can- 
didates in the area selected to represent the spatial position 
of the leading arm in K/M-giants were overdense in velocity 
space. They find a similar trend in velocity space, but with far 
more outliers, probably due to the higher level of contamina- 
tion in the BHB star sample and the broader selection box for 
the leading arm. From our spatially selected sample, we mea- 
sure an average velocity dispersion of Sgr stars of 37kms"'. 
We further select stars within 1 cr of this velocity overdensity 
to be in our 'clean' sample of ^ 70 Sgr BHB stars; such a sam- 
ple suffers from the inhomogeneous angle coverage and dis- 
tance limitations inherent to spectroscopically-selected SDSS 
BHB stars, but has the advantage of high fidelity. 

Using the spatially selected Sgr BHB star sample we exam- 
ine the observed bifurcation of the Sgr stream on the sky and 
its implications for the distances of the two branches. Dif- 
ferent distances of the two branches of th e stream have been 
rep orted by e.g., Belokurov et al. (2006ai. lWnny et al.l (120091) 
andlNiederste-Ostho lt et al.i (i2010i) and are also predicted by 
models (Penarrubia et al. f l2010bl) . This bifurcation is of par- 
ticular interest for the understanding of the formation history 
of the Sgr stream. It has been a challenge for models to repro- 
duce and explain the origin of this feature. It was proposed to 
have originated from wraps of different age (iFellhauer et alJ 
120061) . but also intrinsic propertie s of the progenitor were 
offered as explanations. Recently, Penarru bia et al.l (l2010bh 
presented a model in which this bifurcation is reproduced 
when a rotating disk galaxy is assumed as the progenitor 
There is no indication of a spatial bifurcation on the sky in 
our BHB star sample therefore we use MSTO stars to de- 
fine the two branches on the sky. We measure the mean 
and width of the two parts along the leading arm and find 
only a small 1-2 kpc offset between the two branches over a 
wide orbital angle range and no clear systematic separation in 
the se nse that one is always clearly closer in th an the other 
(cT^^ lYanny etani2009t ICorrenti et al.ll20Tol) . lYanny et al.l 
(120091) reported a trend by visual impression that the branch 



at higher galactic latitude tends to be a t closer distances than 
the branch at lower galactic latitude. ICorrenti et al.l (120101) 
found similar to slightly higher offsets between the distances 
of the two branches compared to this study, but here again 
no clear trend is seen for one branch being always at closer 
distances than the other Although the actual values for mean 
distances of the two branches in our data set seem to be quite 
sensitive to the precise location of the separation cut on the 
sky, we clearly do not see a se paration on the level shown in 
iNiederste-Ostholt et al.l (120101) for the same range of positions 
along t he stream or predicted by the models (Penarrubia et al] 
l2010lj) . It is worth noting that this distance separation ap- 
pears here to be much smaller than the separation of the two 
branches on the sky which is at a 10° level. This strong dis- 
crepancy between the small line-of-sight separation and the 
much larger separation perpendicular to that might be chal- 
lenging to reproduce in the models. 
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